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Abstract 



This is a very elementary introduction to the Heisenberg (XXX) quantum spin 
chain, the Yang-Baxter equation, and the algebraic Bethe Ansatz. 



1 Introduction 



Like the quantum harmonic oscillator, the Heisenberg quantum spin chain is one of the 
fundamental models of physics. Indeed, it is simply formulated, it describes real systems 
which are experimentally accessible [], it has a rich and elegant mathematical structure, 
it is the prototype of all integrable models, and it has close connection to integrable and 
conformal field theory. Here we review the formulation of the model, and we describe the 
elements of the so-called algebraic Bethe Ansatz approach for diagonalizing the Hamiltonian. 
The aim of these notes is to enable beginning graduate students to start exploring the vast 
literature on this subject. 



2 Heisenberg spin chain 

We formulate the spin 1/2 isotropic Heisenberg (XXX) quantum spin chain first for one site, 
then for two sites, and then finally for N sites. 



2.1 One site 

Actually, the one-site problem is too trivial: one cannot write down the Heisenberg spin 
Hamiltonian for only one site. Nevertheless, it is important to identify the set of observables 
and the vector space on which these operators act. 

For the one-site problem, a basis for the observables consists of the identity matrix I and 
the familiar Pauli spin matrices a, that is 

i=i : ; i , *■=! ; : i , °* = (° ~* i , ° z = \ : \ i • w 

These operators act on a two-dimensional complex vector space V = C 2 , whose elements we 

represent by two-component vectors x = I 1 ] , with Xi G C. 

\ x 2 

2.2 Two sites 

We generalize to the case of more than one site using the notion of tensor product. For the 
2x2 matrices ^4 = ( 11 12 ] and B = ( 11 12 ] , the tensor product A®B is defined 
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1 See, e.g., F. Essler's contribution to these Proceedings. 
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For the two-site problem, the basic observables are the spin operators at each site, <7i 
a ® I and a 2 = I <E> a. That is, Q 
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2 I like to draw the horizontal and vertical lines to help organize the matrix elements; they have no other 
significance, and can be omitted. 

3 For visual clarity, I write only the nonzero matrix elements in the tensor product. 
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Note that 



(5) 



for any i,j G {x,y,z}. 

These operators act on the tensor product space V <S> V, with elements 
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This is a good place to introduce the very important permutation matrix V, which is 
defined by 
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In view of Eq. (H) , 
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and so indeed V permutes the factors in the tensor product. Evidently, V 2 = I ® I. Note 
also that 



V fi\ V = a 2 , V 02 V = fi\ . 
The two-site Heisenberg spin Hamiltonian is 4 

#12 = J -{ox-02-\®V)= J - (pfo% + a\a y 2 + ala z 2 - I® I) 



(9) 



We subtract a constant (proportional to the identity matrix) for later convenience. 
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= J - (a x ® a x + a y <g> ^ + a z <g> a 2 - I ® I) 
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We observe that 



H 12 = ^(V-I®I) . 



(11) 



That is, apart from multiplicative and additive constants, the two-site Hamiltonian is given 
by the permutation matrix. 

We consider now the problem of diagonalizing the Hamiltonian, 
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One can verify by inspection that the solution is given by 
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(12) 



(13) 



The first three states have the same energy E — 0, and the last state has energy E = — J. 
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This result can be easily understood from the su(2) symmetry of the model. Indeed, the 
total spin operators 

S = ^(a®I + I® a) (14) 

are generators of a reducible 4-dimensional representation of su(2). That is, there exists a 
unitary matrix U such that 



S( S =i) 






S(s=o) 



(15) 



where S(s=i) and S(s=o) generate irreducible representations of dimension 3 and 1, respec- 
tively. Moreover, 

S 2 = \(a®I + I®a) 2 
4 

= \ (a 2 <g> I + 2B ® a + I <g> a 2 ) 
1 3 

= -a®a + -1(8)1, (16) 
since a 2 = 31. Hence the two-site Hamiltonian can be expressed in terms of S 2 , 

H 12 = a-l®l) = ^{S 2 -2I(8)l) . (17) 

Recalling the well-known fact S 2 \S,S Z ) = S(S + 1)\S,S Z ), we see that 

H 12 \S,S Z ) = ^{S(S + 1)-2)\S,S Z ), S Z = -S,...,S; S = 0,1. (18) 

In particular, for S = 1 the energy is E — 0, and for 5 = the energy is E = — J, which is 
the result we found earlier. 

Note that for J > (antiferromagnetic), the ground state is a spin singlet (S = 0) state; 
while for J < (ferromagnetic), there is a degenerate ground state. 



2.3 N sites 

For the iV-site problem, the basic observables are a n , n = 1 , 2 , . . . , N, defined by 

1 n N 

I J- I 

O n =\ (g) • • • ® I® <7 ®I ® • • • ® I . (19) 



These are operators on 



1 n N 

V ®---®V ®---®V (20) 
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which act nontrivially on the n space, and trivially on the rest. 

There are two possible topologies for a one-dimensional chain: open or closed. Corre- 
spondingly, there are two Heisenberg spin Hamiltonians describing nearest-neighbor interac- 
tions: 

N-l 

H=J2 H n,n+i (open) (21) 

71=1 

or 

N-l 

H=Y^ H n,n+i + H NA (closed) (22) 

n=l 

where the two-site Hamiltonian is 

Ha = jfa- % - 1® N ) ■ (23) 

For simplicity, we shall henceforth focus on the closed chain. 

How can one diagonalize the Hamiltonian? Using a computer is not practical for large 
values of N. Indeed, H is a 2 N x 2^ matrix. Moreover, direct numerical diagonalization will 
give all eigenstates, while often one is interested in only low-lying states. 

A very elegant alternative approach which overcomes (at least in part) both of these 
difficulties is the Bethe Ansatz, of which there are several variants. We focus here on the 
algebraic Bethe Ansatz. 5 An essential element of this approach is a matrix R(X) that is 
a solution of the Yang-Baxter equation, which we discuss in Section 3. As we shall see in 
Section 4, the two-site Hamiltonian ( f23|) can be expressed in terms of the R matrix. The 
fact that the R matrix satisfies the Yang-Baxter equation leads to the "integrability"of the 
model, i.e., that the model can be solved by Bethe Ansatz. 



3 Yang- Baxter equation 

We consider the R matrix Q 

R(X) = XI®I + iV 

5 The method pioneered by Bethe Q is now known as coordinate Bethe Ansatz. The algebraic Bethe 
Ansatz was developed in St. Petersburg Q - ||. It is perhaps not as powerful as coordinate Bethe Ansatz, 
but it is more transparent. Other related approaches include analytic Bethe Ansatz W and Baxter's Q- 
operator method ||. A rather different approach, which is formulated directly in the thermodynamic limit 
and avoids Bethe Ansatz, is being developed in Kyoto ||. 

6 The notation 0] which we are using, while widely used, has not been universally adopted. In particular, 
some authors (e.g., ||, §, |J) call R{\) what we call R(X) = VR(X). Given an R matrix, one can easily 
determine which of the two conventions is being used by evaluating it at A = 0: for us R(0) ~ V, while for 
them R(0) ~ I <g> I. Of course, their Yang-Baxter equation looks somewhat different from ours. 



6 



/ X + i 



X 



X 



X + i 



where 



a = X + i , b = X , c = i . 
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We regard -R(A) as an operator which acts on V <8> V. The variable A is called the spectral 
parameter. 

We wish to show that this matrix is a solution to the so-called Yang-Baxter equation. 7 
To this end, we define i?i2(A) by 
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V 

Evidently, this is an operator on V <8> V <S> V, which acts nontrivially on the first and second 
spaces, and trivially on the third. 

Similarly, we define #23 (A) by 
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7 For the case of the simple R matrix (pi|), the approach we shall follow is unnecessarily tedious - there 
are much faster ways to show that this matrix is a solution of the Yang-Baxter Eq. (|3l|). However, our 
approach has the benefit of being explicit and straightforward, and it works for any case. 
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V 

This operator acts nontrivially on the second and third spaces, and trivially on the first. 

Finally, we wish to define #13 (A), which acts nontrivially on the first and third spaces, 
and trivially on the second. This requires only slightly more effort: 

#13 (A) =V 23 #12 (A) V 23 , (28) 

where #12 (A) is given by Eq. fl26|) , and V23 is given by 

(l \ 
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One can now easily obtain 
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The Yang-Baxter equation is 

R 12 (X - A') R 13 (X) i? 23 (A') = R 23 (X') R 13 (X) R 12 (X - A') . (31) 

It is now a tedious but straightforward exercise in matrix multiplication and algebra to show 
that this relation is satisfied. With the help of a symbolic manipulation program such as 
Mathematica, the task can be accomplished in minutes. f\ 



4 Transfer matrix 

In this Section, we show that the two-site Heisenberg Hamiltonian can be expressed in terms 
of the R matrix, and that the model is integrable by virtue of the fact that this matrix is a 
solution of the Yang-Baxter equation. However, it is more convenient to work "backwards" : 
using the R matrix, we first construct the so-called transfer matrix - a one-parameter com- 
mutative family of operators acting on the space of states of the Heisenberg spin chain. We 
then verify that the Heisenberg Hamiltonian is among this family of commuting operators. 

The key step in this program is to introduce the so-called L operators f\ 

Lon(X) = 

z 

( n, R \ 

n = l,2,... ,N, (32) 
which act on so-called auxiliary (0) and quantum (n) spaces. That is, L Qn is an operator on 

1 n N 

V ®V (33) 
which acts nontrivially on the th and n th spaces, and trivially on the rest. Moreover, 

1 n T 

a n = I®-..®I®a®I®-.-®I, 

Pn = I®-"-®I®/9®I®---®I, 

7„ = !®---®I®7®I®---®I, 

6 n = I®---®I®5®I®---®I, (34) 




8 The R matrix ( |24| ) is the simplest solution of the Yang-Baxter equation. Many more solutions are known. 
See, e.g., |, @. 

9 The shift in A is made in order that the Bethe Ansatz equations (see Eq. (^Tj) below) have a symmetric 
form. This shift is not essential. 
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where 




Evidently, a n , /3 n , j n: 5 n are operators on 

1 n N 

V ® • • • ® y ®---®V (36) 

which act nontrivially on the n th space, and trivially on the rest. The vector space fl36|) 
is precisely the space of states of the Heisenberg quantum spin chain (p0|). Note that the 
operators at site n commute with the operators at any other site n' ^ n. 

The algebra of the operators a n , (3 n , j n , 5 n is encoded in the relation 

-Roo'(A — A ) Lon(^) L / n (\') = L ' n (A ) Lon(A) i?oo'(A — A') , (37) 

which immediately follows from the Yang-Baxter Eq. ([31]) upon replacing 1 —>■ 0, 2 — > 0', 
and 3 — > n, and also making the shifts A — ► A — | and A' — > A' — \. 

The monodromy matrix T (A) is defined as the following product of L operators 

T (A) = L 0N (X)---L 01 (X) 

f a N f3 N \ (a. x fa 
\ In $n J \ 7i <5i 

The monodromy matrix obeys the fundamental relation 

R 00 ,(\ - A') T (A) T ,(\') = T ,(A') T (A) i? 00 ,(A - A') . (39) 

We prove this result for N = 2. The left hand side (LHS) is then 

LHS = R 00 ,(\ - A') L 02 (A) Lqi(A) L ^(A') L a(A') 

= -Roo'(A — A') L 02 (X) L i 2 (\') Lqi(X) L 'i(A') 

= L '2(A') L 2(A) -Roo'(A — A') ^oi(A) Lo'i(A') 

= L '2(A ) L 2(A) Lo'i(A') L i(\) i?oo'(A — A') 

= L Q , 2 (A') Loa(A') L 02 (A) L 01 (A) ^(A - A') = RHS . (40) 

10 As is customary, we often suppress the quantum-space subscripts. 
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In passing to the second line, we have used the fact that L 0n commutes with Low f° r n ^ n', 
and in passing to the third line we have used the fact (j37D - This proof immediately extends 
to the case of general N. 

The transfer matrix t(X) is defined by tracing over the auxiliary space 

t(A)=tr T (A). (41) 

It acts only on the quantum spaces fl3"6"|). The transfer matrix constitutes a one-parameter 
family of commuting operators 

[f(A),t(A0] = O. (42) 

Indeed, multiplying the fundamental relation (P5|) on the right by the inverse R matrix, and 
taking the trace over both auxiliary spaces and 0', we obtain 

tr«y i?oo'(A - A') T (A) T ,(X') R 00 ,(X - A')" 1 = tr 00 , T ,(X') T (A) . (43) 

Using the cyclic property of the trace, we obtain 

trw T (A) T ,(X') = tr 00 < T ,(A') T (A) , (44) 

which implies the result (fl2|). 

The transfer matrix contains the momentum operator P of the (closed) quantum spin 
chain. Specifically, 

p = iio g r N t( l ~). (45) 

We show this for the case N = 2. We first observe 

t( 1 -) = tioToC-) = tr i? 02 (0) Rox(0) = t 2 tr V 02 V i = ~Vu , (46) 

where we have used the fact R(0) = iV (see Eq. (|24])) and the identity XxqAqVqu = A n . It 
follows that 

t(^) X x tC-y 1 = V12 X, V 12 = X 2 , (47) 

where X n is any operator at site n. One can show for general N in similar fashion that t(|) 
is the one-site shift operator, 

t(^)X n t(^)- 1 = X„ +1 , n = l,...,N. (48) 
We recall from elementary quantum mechanics that 

e iaP x e -iaP =x + a, (49) 
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where P and X are the momentum and position operators, and a is a c-number. By analogy, 
we define the momentum operator for a spin chain by e lP = i~ N t(^), which implies ([45]) . 

The transfer matrix also contains a closed-chain Hamiltonian 

H= 2VdX l0gt(A) a=- ~ Nim ) = S Hn > n+l +Hn ^ (50) 



71=1 

with the two-site Hamiltonian 

^ = ^(^•^•(0)-^) , (51) 

where the prime denotes differentiation with respect to A. This is easily verified for N = 2 
using (f46[ ) and also 

t'C-) =tr o i4(0) i2oi(0)+tr o i2o 2 (0) ^(0) = 2iR' 12 (0) . (52) 



Moreover, keeping in mind the result ( |lT|) and the expression ( p4|) for our i? matrix, we 
see that the two-site Hamiltonian ([51]) coincides with that of the Heisenberg chain (p3|) . In 
short, the Hamiltonian (|50|) is precisely the Heisenberg Hamiltonian (p2|). 



Combining the results (|42j) and (0), we see that the Hamiltonian commutes with the 
transfer matrix, 

[H,t(X)) = 0. (53) 
This relation together with ([42"D signal the "integrability" of the model. 



5 Algebraic Bethe Ansatz 

We finally turn to the problem of diagonalizing the Hamiltonian. In fact, we shall diagonalize 
the transfer matrix. The strategy is to identify certain creation operators with which to 
construct states, as one does for the simple harmonic oscillator. These operators, however, 
depend on parameters; the states are eigenstates of the transfer matrix if the parameters are 
solutions of a set of equations, namely the Bethe Ansatz equations. 

The monodromy matrix T (A) is a 2 x 2 matrix in the auxiliary space, whose elements 
are operators on the quantum space V® N which we denote as follows 

1 C(A) D(X) 1 
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A set of algebraic relations among these four operators is encoded in the fundamental relation 
fl39f). The relations include 

[B(X),B(X f )]=0, 

D(\)B(\>) = ^^B(X')D(X)- C ^-^B(X)D(X'), (55) 



where a , b , c are given by Eq. (^5|) . 

Let uj + be the ferromagnetic vacuum state with all spins up, 



N 

which is an eigenstate of ^4(A) and D(X), and which is annihilated by C(X), 

A(X)u + =(x+-\ lu + , D(X) cu + =(x- % -\ lu + , C(X)u + = 0. (57) 

We use the operators B(X) as creation operators to construct the so-called Bethe state 

|A X) ... ,X M ) = B{X 1 )---B{X M )u J+ . (58) 

Using the algebraic relations fl55| ) and the properties flo"7| ) of the ferromagnetic vacuum state, 
it can be shown that the Bethe state is an eigenstate of the transfer matrix t(X) = A(X) + 
D(X), 

t(X) |Ai , . . . , A M ) = A(A ; Ai , . . . , X M ) |A X , . . . , A M ) (59) 
with the eigenvalue 

/ • \ N M \ \ ■ / ■ \ N M \ \ i 
A(A;V...,A„H A + i n^I^H II (») 

x 7 a=l x 7 a=l 

if {Ai , . . . , Xm } are distinct and obey the Bethe Ansatz equations 

K±,\ =U K^ ±1 a=1 ,.., M . (61) 

The proof is immediate for N = 1, and is not much more difficult for general N. It is 
well- described in many papers (see, e.g., @,@,@), so we shall not repeat it here. 
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In particular, one can now obtain the momentum and energy eigenvalues using Eqs. (|45|), 
and ©, 

i M /i ±«\ j M i 

An important feature of the Heisenberg spin chain is its su{2) symmetry. One can show 
H that the transfer matrix commutes with the total spin operators, 

1 - 

= 0, S =-J2°n. (63) 

n=l 

Moreover, the Bethe states (|58|) are su(2) highest weight states 

S + |Ai,... ,A M ) = 0, S ± = S x ±iS y , (64) 

and are eigenstates of S z with eigenvalue 

S z = ^-M. (65) 

Since S = S z > 0, it follows that M < N/2. Lower weight states are obtained by acting on 
the Bethe states with the lowering operator S~ . 



S,t{\) 



6 Now the hard work begins 

The S = 1/2 isotropic Heisenberg spin chain which we have discussed is the simplest inte- 
grable magnetic chain. Many generalizations are possible. For instance, using appropriate 
R matrices, one can construct anisotropic (XXZ and XYZ) chains, chains with spins in 
higher-dimensional representations (S = 1 , 3/2 , . . .) or with combinations of different spins, 
and chains with spins in representations of higher-rank algebras (su(A/") ,...). Using also 
so-called K matrices, one can construct integrable open spin chains [JTTj] . Bethe Ansatz 
equations have been obtained for many such models. 

Obtaining Bethe Ansatz equations (BAE) for an integrable model is in a sense "half" the 
solution of the model - the easier half at that! The difficulty is that quantities of physical 
interest (e.g., energy) are expressed in terms of solutions of BAE; and unfortunately, BAE 
are generally hard to solve. Various clever methods have been developed for extracting 
from these equations information about the thermodynamic (N — > oo) limit. Invariably 
the methods involve some assumptions about the nature of the solutions, such as the string 
11 Numerous additional references can be found in flql. 
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hypothesis (see, e.g., ||). The types of results which have been obtained include spectrum 
of low-lying states, scattering matrices, finite-size corrections, and thermodynamics (finite 
temperature and magnetic field). Progress has been made on the computation of correlation 
functions. See, e.g., and other contributions to these Proceedings. 
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